Simple theory of the measured current through quantum dots 
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A simple theory of the detected current I{t) flowing through charge qubits — quantum dots — 
is proposed in terms of standard continuous measurement theory. Applied to a double dot, our 
formalism easily confirms previous results on quantum Zeno effect, driven by growing ammeter 
performance 7. Due to the transparent formalism, we can calculate the exact fluctuation spectrum 
S{uj) of the detected current, containing a significant Lorentzian peak near the Rabi frequency of 
the double dot. 
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Back-action of measurement on the measured object is 
an emblematic property of quantum systems. Quantum 
Zeno effect (QZE) - suppression of the object's coher- 
ent internal dynamics by measurement - is one of the 
marked scenarios for that back-action. In this Letter we 
are going to discuss the effect in the framework of un- 
sharp measurements, for which there is a well-developed 
theoretical framework called time-continuous weak mea- 
surement theory (cf. Refs. 12] and references therein). 
The application we have in mind is the case of a double 
quantum dot (DQD), which is a semiconducting nanos- 
tructure, available in high quality due to massive progress 
in experimental technology. 

There is a growing culture of indirect measurements 
on nanostructures by means of Coulomb-coupled quan- 
tunipoint contacts, single-electron transistors, or DQD's 
i-i. QZE is one of the effects studied from the begin- 
ning pI (cf. also 0), and the concept of time-continuous 
measurement has penetrated the field ^, ^ for a long 
time. All those studies assume the - sharp or unsharp - 
detection of the number N{t) of electrons that have tun- 
neled through the nanostructure, and the current I(t) is 
defined as the stochastic mean {dN{t)/dt). The present 
work differs from those studies in assuming that detec- 
tion is done by a tool usually considered as fully classical: 
an ammeter of high performance, monitoring the time- 
dependent current I{t) flowing through that device. 

The main parameter of the theory is detection perfor- 
mance (or detection strength), defined as 



7 = (A0-i(A/)- 



(1) 



where At is the time-resolution (or, equivalently, the in- 
verse bandwidth) of the ammeter and A/ is the statisti- 
cal error characterizing unsharp detection of the average 
current in the period At. The accuracies of commercial 
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ammeters reach pA « 10^ electron/ s, at a bandwidth of 
10** Hz. Then 7 — 10"^" s, which should be compared to 
the time scale of internal coherent dynamics of a DQD, 
characterized by the Rabi frequency Vl « 10"'^'' Hz, which 
is also the order of magnitude of the steady-state current 
measured in electron/ s units, viz. I /e. After all, ~ 1 
can be reached through standard instrumentation. Be- 
low we are going to show that this suffices for observing 
a continuous-measurement version of QZE. 

We rely upon the standard Markovian approximation, 
tracing out environmental variables, and treating trans- 
port in terms of the reduced density matrix p of the DQD 
and a corresponding effective current operator /. Then 
the stochastic mean of the detected current is obtained 
as the quantum mechanical subsystem average Q 



(2) 



Continuous measurement theory starts from this point, 
and accounts for a double action of the ammeter of finite 
performance: 

1. The quantum mechanical back-action on the con- 
tinuously measured quantum system, our main concern 
here, induces loss of coherence between eigenstates be- 
longing to different eigenvalues of the measured quantity; 
the rate of decoherence being proportional to the detec- 
tion performance 7. Under well-defined conditions [11, Q 
it can be represented by adding a Lindblad-type term to 
the master equation: 



dp 
'dt 



cp~l[i, [i,p]]^c'p, 



(3) 



Our analysis is based on the solution of this extended 
master equation, in which the Lindblad supermatrix C 
provides a Markovian description of decoherence and 
damping, as usual. 

2. Along with the quantum effect we expect to observe, 
inaccurate measurement generates fully classical white 
noise in the measured current, superposed on the average 



I{t) ~ {I(t)) - gw{t). 



(4) 
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Here w{t) is white noise on the scale of the Marko- 
vian dynamics (i.e., a fiat spectrum noise of band- width 
much larger than the typical frequency range correctly 
accounted for by the Markovian approximation). The 
amplitude of the noise is related to the detection perfor- 
mance 7 according to an inverse square-root law: 



g = i/%/7- 



(5) 



Before doing the calculation, we observe that a genuine 
ammeter, unlike an electron counter, smears out dynam- 
ical details of shot noise originating from incoherent tun- 
neling of individual electrons through the external bor- 
ders of the DQD, while retaining time structure due to 
internal coherent dynamics (Rabi oscillations). Accord- 
ingly, we omit shot noise from the present treatment; 
hence the calculus becomes surprizingly simple and ro- 
bust, as can be easily checked for a single dot where the 
power spectrum is exactly known from the iV-resolved 
calculations p^ . 

Let us introduce the stationary solution poo of the ex- 
tended master equation ([3]): 



C'poc = , 

as well as the stationary current: 

loo = (/(i))oo = (/)p. 



(6) 



(7) 



where (.)oo stands, in general, for stochastic mean of cur- 
rents detected on p^o ■ We also define the Heisenberg op- 
erator of the current for < > 0: 



dijt) 
dt 



{c')n{t) , /(o) 



(8) 



In continuous measurement theory, the stationary corre- 
lation function of the fluctuating detected current turns 
out to be: 



mm)oo {i{t))r_ 

= Re{i{\t\)i)p^ 



>{m)c 



(9) 



The nonsingular term on the r.h.s. follows from the stan- 
dard expression (/(t)/(0))oo = R'e(/(|t|)/)p^ valid for 
bulk quantum systems where the detector noise can be 
neglected. For our nanostructure, the singular term at 
t = comes from the white-noise fluctuations, originat- 
ing from unsharp measurement, i.e., detection noise, in- 
troduced in eq. ^ . We have thus expressed the station- 
ary correlation function of the detected classical current 
I{t) in terms of the quantum correlation of the Heisen- 
berg current I{t). The spectral density of the detected 
fluctuations will be defined as the Fourier transform [§| 
of the correlations resulting in the form: 



Sico) = / (/(t)/(0)), 



27tIIS{lu) , (10) 



which for high frequencies approaches the measurement- 
added white noise value 



For very high values of ammeter performance 7 that noise 
level may get below the "shot noise limit" /oo/2. In that 
case ammeter and charge counter measurements may dif- 
fer in characteristic ways not controlled by the present 
method of calculation; however, this is not our concern 
here, since the effect we envisage remains in the well- 
treated range of moderately high 7. 

Now we turn to the specific features of the double 
quantum dot, consisting of two potential wells ("dots": 
left and right), with an internal barrier allowing coher- 
ent tunneling between them, and two external barriers 
allowing incoherent tunneling between each dot and its 
joining lead. The Markovian approximation used is based 
on the assumption thermalization in leads L and R being 
the fastest dynamical process present. 

Because of intra-dot Coulomb blocade, at low temper- 
ature in each of the dots there can be but one electron 
in the ground state, or none. Strong inter-dot Coulomb 
repulsion further reduces the set of available orthogonal 
basis states to the following three: 

|0) = |0,0), |i) = |l,0), |i?) = |0,l). (12) 

On the above basis, we introduce the following absorption 
and emission operators: 

aL = |0)(i|, dR^\0){R\, (13) 

as well as the charge operators: 

hL = \L){L\, hR^\R){R\, fi^nL + nR = l- 

With the above definitions, the Hamiltonian is this: 

^0 = h^fiL - nn) + n{d{dR + H.C.) . (15) 

Neither the Hamiltonian, nor any of the relevant ob- 
servables (in particular, the current - see below) have 
off-diagonal elements connecting the state |0) to the rest 
of the reduced Hilbert space of the DQD, imposing an ef- 
fective charge superselection. Therefore the off-diagonal 
elements OL, Oi?, LO, i?0 of the matrix p should be set to 
zero identically: 



(14) 



poo 
Pll Plr 
Prl Pbr 



(16) 



5((X3) = 5' = 1/7- 



(11) 



The Markovian approximation yields two irreversible 
processes that modify the intrinsic Hamiltonian dynam- 
ics of the dots: tunneling of an electron from lead L to 
the left dot at rate F^, and tunneling of an electron 
from the right dot to lead R at rate Vr. The rates 
Tl , Ffl depend on the details of the total Hamiltonian 
dynamics of L-|-R-|-dots, and we shall take their value for 
granted. The two processes contribute to two currents, 
II — Fi(l — h) and Ir = Trjir. Taking into account 
the Ramo-Shockley effect of fast screening in the leads 
[Tl| . on the slow Markovian time scale there is a single 
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time-local operator of the current flowing from L through 
dots and R to the ammeter: 

i.£^(1^4^±I-^.if'o'o . (17) 

2 ^ \ Tn J 

The Lindblad supermatrix C appearing in eq. ([3]) con- 
tains the modification of Hamiltonian dynamics due to 
the two tunneling currents Il/r- It is easy to work out: 

Cp=-i[Ho,p]+ (18) 
+^L{alpaL ~ HaLO-l, p}) + Ti^XaHpa^ - HaflOfl- P}) 

where Lamb shifts have been included in the energy split 
S of Hq. Note that the continuity equation is always 
satisfied [ll,[i|: 

iR-iL+C^n = , (19) 

where £^ is the adjoint of C. 

Quantum measurement back-action enters through the 
modified Lindblad supermatrix C introduced in eq. ([3]). 
Then the stationary state is obtained by solving eq. ([6]), 
in which eqs. p^ . p?]) and are used as input. 
The calculation is standard; surprisingly, we obtain ex- 
actly the same equations as those derived by Gurvitz [3] 
for the closely related but still different model of a DQD 
observed by a point-contact charge detector, if Gurvitz's 
clicking rate Td is replaced by our combination ^Tj^/W. 
We think the reason is that as long as Markovian ap- 
proximation is valid, the effective current operator / is 
pinned down to the single-dot occupation operators ul/r 
through eq. P7)) . Then a point-contact charge counting 
device, like the one discussed by Gurvitz, and a commer- 
cial ammeter, as we suggest, have no more freedom than 
to couple to the DQD in the same way. 

The stationary electron number current reads: 

J ^ TLTRff'jl + y) 

°° S^Tl + {2rL + TR)n^l + y) + ^TlTKI + y f ' 

(20) 

where y = 7ri^/16. The above formula contains the way 
the Quantum Zeno Effect appears with growing measure- 
ment performance 7 (see Fig. 1): in all cases, /qo ^ 
for 7 — !■ cx). By inspection of the density matrix we 
learn that this is due to damping of the coherent interdot 
transport, causing increased occupation of the left dot, 
blocking the way of new electrons to enter. With no bias, 
5 — the reduction of current is monotonous; with suf- 
ficiently strong bias, 5/Q, > 2, however, there is a range 
of small performances where the current increases with 
growing 7 ("anti-Zeno effect"). Apparently, asymmet- 
ric occupation induced by decohering measurement and 
that caused by bias are competing with each other. The 
whole effect is small though: for <5 = 0, as the measure- 
ment performance passes the shot- noise limit 7 « 2/ loo, 
it reaches a 3% reduction of current. 




FIG. 1: Measured stationary current through an unbiased 
double quantum dot, as a function of ammeter performance, 
for symmetric tunneling rates Tl ~ F^ — Q (eq. ([l])) (solid 
line), and asymmetric ones with Fr = 5F_l ~ ^ (dashed line). 
The shot-noise limit is at ^Q, = 6.7. 

Of particular interest for the experiment is the case of 
asymmetric tunneling rates, Tr > Tl, resulting in much 
stronger QZE, reaching 30% for the example displayed 
in the figure. In this case the right dot is depleted with 
respect to the left one, which causes higher sensitivity to 
measurement-induced blocking of the left-to-right tun- 
neling. Asymmetry of the opposite sign has no marked 
effect on QZE. 

Our simple theory presents the best of its performance 
in evaluating unsharp measurement back-action effects 
on the detected current fluctuation spectrum. Equations 
([9|) and (fTO| set the framework for the calculation; in 
particular, it is the operator correlation function appear- 
ing in eq. ([9]) that has to be evaluated by means of our 
Markovian master equation. We outline the unbiased 
case S = 0; then 

C'i\L)(R\ + \R){L\) = -^,TRil+^,TR/16)i\L)(R\ + \R){L\) 

(21) 

hence the real part of the LR component of the Heisen- 
berg current I{t) defined by eq. ([8]) remains zero. So the 
matrix I(t) has four independent components instead of 
five. We can thus replace eq. © by di{t)/dt = C\H{t) 
where £j is the restriction of C on the four-dimensional 
subspace explored by I{t). Note, furthermore, that C is 
degenerate also because C poo = and this degeneracy is 
inherited by C'j, too. Therefore the characteristic equa- 
tion |£j -t- A| =0 has a trivial zero root, and the rest of 
the characteristic equation is only cubic. For the special 
case Fi = F/f = r2 it turns out to be: 

A^- P + 2^) i7A^+ fe + n n^x^ ('1 + 2^^^ ^0 . 

V2 32 y V 16/ V 2 32/ 

(22) 

It has one real positive root A = Aq and two complex 
conjugate roots A = Ai ± iuR with positive Ai, imply- 
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FIG. 2: Noise spectrum of the measured double-dot current, 
as a function of the ammeter performance 7. 

ing that the solution I(t) consists of an exponentially de- 
caying and an exponentially damped oscillatory part, the 
latter oscillating at the damped Rabi frequency cor « 2fl. 
Let us have a look at the structure of the corresponding 
fluctuation spectrum PH)) : 

"^"^'-Lj^ + A2"^(l^ - ivjiy + + + A? ■ 

(23) 

One can obtain closed expressions for the parameters of 
the above expression. We present the results for low val- 
ues of the dimensionless combination ^fl: 

ujr/2Q = 1.025 + 0.0027O-HO(7^^)^ 
Ai/0 =0.52 + 0.0U"/n + O{'ynf. (24) 



This spectrum has one local minimum at a; = and 
two wide peaks at the Rabi frequency iw^j, the widths 
of the peaks are cca. one fourth of the Rabi frequency, 
and on increasing the performance of the ammeter their 
amplitudes start to decrease. For high values of 7 the 
Rabi peaks get overdamped and merge into a single peak 
around uj = 0, see Fig. 2, where background noise has 
not been included. We stress that to observe the peaks 
themselves a good ammeter (high 7) is needed, since a 
poor ammeter introduces too much background white- 
noise g"^ = 1/7 to see anything of the coherent peak 
structure. 



As a conclusion, we have derived explicit expressions 
for the back-action of an ammeter on measurable char- 
acteristics of a double quantum dot. The back-action 
has the character of quantum Zeno effect, counteract- 
ing coherent internal motion of the object, revealed both 
as reducing the mean transmitted current, and damp- 
ing - eventually, overdamping - Rabi oscillations, as ob- 
served in the noise spectrum. In the system studied, mea- 
surement back-action takes place through decoherence- 
controUed asymmetry in the occupation probabilities of 
the two coupled quantum dots. Other effects influencing 
the same asymmetry may compete with the effect, even- 
tually producing ranges of anti-Zeno effect. Appropri- 
ately tuned asymmetric tunneling barriers may strongly 
enhance the possibilities of observing the effect. 
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